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The Elliptic Inequalities in the Lunar Theory. 

By Ernest W. Brown, Fellow of Christ's College, Cambridge. 



The object of this paper is to give a general solution in series of the 
problem "A system of three bodies is in motion in one plane, the first is revolving 
about the second and is disturbed from its elliptic orbit by the third. The 
third body is supposed to be of infinite mass and to be moving in a circle of 
infinite radius with a finite angular velocity. Given the relative positions of the 
three bodies at any one time, to find their relative positions at any other time." 
The differential equations which express the motion of the first body relative to 
the second have been given by Dr. Hill in Vol. I of this journal. The second 
body is taken as the origin, and the motion is referred to rectangular axes mov- 
ing with the angular velocity of the third body which always lies on the x-axis. 
The equations are 

*-*£+(£-*■>=•■ 

d?y , „ , dx , u 

df+^dt +->» = ' 

where y. is the sum of the masses of the first and second bodies and n' the angu- 
lar velocity of the third about the origin. 

One integral of these equations is immediately obtained. It is 

(&)'+(#) = £+*«-«'• 

In chapter II of his memoir just referred to, Dr. Hill transforms these 
equations to the imaginary variables x ± yV — i, and by the use of this, the 
Jacobian Integral, brings them to two simultaneous differential equations of the 
second order and second degree. 
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A particular solution is then found which can be expressed by an infinite 
trigonometrical series with a single period, giving as an orbit relative to the 
moving axes, a closed symmetrical curve which, for small values of the ratio of 
the mean motions of the first and third body, is an oval. I have given the 
extension of this particular solution to the case where the third body is at a 
comparatively great but not infinite distance.* The orbit is still a closed curve, 
but it has lost its symmetry about the axis of x. This extension presented no 
special difficulties beyond the question of representing the series obtained in the 
most convergent form. 

When an attempt is made to obtain the general solution, it appears below 
that it will be necessary to represent it by an infinite series with two periods ; 
that is, by a doubly infinite trigonometrical series. The orbit is now no longer a 
closed curve, but it is one which re-enters its orbit after an infinite time, and the 
maxima and minima values of the coordinates lie between determinate limits. 
The chief difficulty lies in the fact that the equations of condition between the 
coefficients of the terms in the series require a relation between the two periods ; 
it is the finding of this relation that entails the trouble. 

Again, since the solution is expressed by an infinite series, it is necessary 
in order that it should have a meaning, that the series should be convergent, 
for this, two of the arbitrary constants must be confined within certain limits. 
We are thereby limited to consider only disturbing bodies whose mean motions 
have ratios to that of the revolving body which do not surpass a certain value, 
and also that the body when undisturbed shall have an eccentricity less than 
some given quantity. It does not seem easy to determine what these limits are 
or whether they exist at all. M. Poincare's remarkable researchesf into the 
methods of the problem of three bodies have to a large extent opened the way 
to the determination of convergency. He has in fact proved that for small 
values of certain of the constants, there are periodic solutions which can be 
represented by convergent series. The convergence of the series obtained 
below will depend in a large measure on the value of c, the ratio of the two 
periods of the general solution. 

The principal part of the ratio c has been shown to be one of the roots of an 
equation expressed as an infinite determinant, and it was this peculiarity that 

* American Journal of Mathematics, Vol. XIV, p. 141. 

t Acta Mathematica, Vol. XIII. Mecanique Celeste, Vol. I. 
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stimulated inquiry into the possibility of the existence of such forms. G. W. Hill, 
Poincare and Helge von Koch have given investigations of their properties. In 
particular, Dr. Hill,* under the assumption that it was possible, actually obtained 
the principal part of c from the infinite determinant which arose in his method ; 
later M. Poincaref gave conditions which, when satisfied, allowed the determinant 
to possess a finite and determinate value, and also proved that it was possible to 
obtain c from it ; and lastly, Dr. von KochJ has shown under what circumstances 
the ordinary transformations used in determinants with a finite number of rows 
and columns can be extended when these latter are infinite in number. 

If we examine the series given by Delaunay in longitude and parallax for 
the coefficients of the several trigonometrical terms, it appears that at any rate 
in the case of the moon, slow convergence, whenever it occurs, proceeds in any 
coefficient along the series arranged in powers of the ratio of the mean motions ; 
also that in terms with actually small coefficients the convergence appears to be 
slower than in the larger coefficients. As pointed out by Dr. Hill, it becomes of 
interest to see how these can be improved either by taking new quantities in 
terms of which to make the expansions or by using numerical values from the 
start. He has examined this question in the case of the principal part of the 
Variational Inequalities, and I have attempted to give convenient expressions for 
the same purpose in that of the Parallactic Inequalities. When the equations 
for the Elliptic Inequalities are examined, it does not appear that there is any 
method of improving the expressions for the coefficients on the same lines 
beyond taking m = n'/(n — n') instead ofn'/n in terms of which to make expan- 
sion. The quantity c, which is itself represented by a slowly converging series 
in powers of ni, is involved in a complex manner, and even if it were not, we 
should simply be brought back to the difficulty of representing c by a quickly 
convergent series. Hence it seems better to use the numerical value of m from 
the outset, keeping all the other constants involved arbitrary ; that is to say, in the 
case treated here the eccentricity is left arbitrary. 

Another consideration justifies this method of proceeding. Of all the con- 
stants used in the determination of the motion of the moon, the ratio of the 
mean motions is most exactly determined, and it is doubtful whether any future 
observations will change it to any appreciable extent. This is not the case with 

•Acta Math., Vol. VIII. t Bulletin de la Soc. Math, de France, Vol. XIV, p. 77. 

I Acta Math., Vol. XVI, p. 217. 
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the other constants. As the slow convergence occurs in the series arranged in 
powers of this ratio, little will be lost by keeping the rest of the constants arbi- 
trary, and the extra trouble on the plan here adopted is very small. If then a 
complete theory of the moon were given on this principle, the coefficient of 
every term would appear as a quickly convergent series arranged in powers of 
the eccentricities, the ratio of the parallaxes, and the sine of the latitude, with 
numerical multipliers only. The coefficients in such a development would have 
a comparatively small number of terms with a high degree of accuracy. 

Another point may also be noted in connection with these methods It 
appears generally that in proceeding to find the coefficients by successive approxi- 
mation, the second and higher approximations are due principally, not to the 
introduction of new terms in the series, but to the second and higher approxima- 
tions of the terms previously used, so that it is not generally necessary to find 
the first approximations of very many of the higher terms in order to get the 
large coefficients accurately. This fact abbreviates the calculations to a large 
extent. 

With respect to the numerical results obtained, no trouble has been spared 
to make them accurate. Equations of verification have been computed for each 
step and also for the final results. The various forms which the equations of 
motion can take enable this to be easily done. In the second portion of the 
paper (to be published hereafter) a value has been obtained for the part of the 
motion of the Lunar Perigee which depends on the square of the eccentricity. 
In order that this value should possess a degree of accuracy which can be relied 
on, it is absolutely necessary that the quantities obtained be verified by inde- 
pendent equations. I have not had the advantage of any assistance in carrying 
out the computations, and have consequently taken special pains to find out the 
best equations for verification, and have computed in certain cases two and even 
three of them to test the same results. 

I. 



The equations under discussion are 




£-*&+(*- 


Sn'^x=0, 


%+«$+*■* 


= 0, 



33 



248 Brown: The Elliptic Inequalities in the Lunar Theory. 

which admit of a general solution when n' is zero and a particular solution when 

n 1 is not zero. In the first case the resulting orbit is an ellipse, and therefore, 

as is well known, x, y can be expressed in terms of the time by means of 

Bessel's functions* 

+ ■» 

x = a\p — J$~ n cos jpg (t — to) , 

+ 00 
00 

1 3 

where — JJ-" = e, b 2 = a 2 (1 — e 2 ) and g is the mean motion. The four 

quantities, a, e, g, t , involve three arbitrary constants, the fourth being used 

in making the major axis of the ellipse coincide with the cc-axis. The particular 

solution of the equations involves only two arbitrary constants, and is of the 

formf 

+ «> 

x — a o /' . a i cos (2* + 1)(« — n')(t — ^i) , 

— oo 

y = a o y< a t sin (2* + l)(n — n')(t — ^), 

— oo 

the arbitraries being n (or a ) and t x (when 4=0,^=1). These two forms are 
distinct unless both e and n' are zero, when they reduce to circular motion. The 
coefficients a t are functions of the ratio n' :n — n'. Since then we have these forms 
for zero values of n' and e respectively, it is assumed that for small values of e 
and n' there is a possible solution of the form 

x = a 2 i 2 p A i<p co S {(2i + l)(n-n>)(t-t 1 )+pg(t-t )\, ,. +oo _x 

y = 0,2,2,4, p sin {(2i+l)(n— n')(t-t l )+ pg(t — < )}, "" 

This solution involves the necessary four arbitrary constants and may be looked 
upon as a general solution of the equations except in so far as e and n'/(n — n') 
are limited in magnitude. This limitation is necessary in order that the series 
may be convergent, but it is to be understood that there is no guarantee that the 

* Tisserand, M6c. Cel. t. I, p. 226. 

tG. W. Hill, Researches in the Lunar Theory, Amer. Jour. Math., Vol. I, p. 130. Unless otherwise 
stated, the references to this memoir are given by the page of Vol. I of this journal. 
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series are convergent for even small values of these constants. At present all 
we know is that the solution with n' zero converges for values of e sufficiently 
small, and the particular solution for sufficiently small values of the ratio 
n' j{n — n'). 

The relation g : n — n' will involve some of the constants. Put then 

g = c(n — »'), 

a=\(2i+l)(t — t 1 )+pc(t — t Q )\. 



Transform the equations of motion by putting u = x + yV — 1 , s = x — yV — 1 , 

and „/ - , ^ n , 

= D, —,=.m; 



V — 1 d __ n n' 

n — n' ' dt 



they become, with the help of the Jacobian Integral (Vol. I, p. 132), 
D 2 (t«) — Du.Ds — 2m (uDs — sDu) + — m 2 (u + sf = G, 
D (uDs — sDu — 2mus) + -— m 2 (u 2 — s 2 ) = . 



We have D (aV- 1) = 2i + 1 +cp, 

and since the coefficients of the various powers of e* v ~ 1 for all values of i, p 
(except the continuation i = p — 0) will be equated to zero, the constants t and t x 
will not enter into the equations of condition thus found. They may be therefore 
omitted in what follows, care being taken that they be inserted when the values 
of the coefficients and of c have been found and the series written out in full. 
Hence we can put 

a = (2i + 1 + cp)(n — n')r, 

where dr = dt. Also write 

so that n v d 

*dj' 

and the assumed solutions become 

u-= 2^, {*+»+•, i 
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Substitute these in the equations (1). For this purpose 

Du = 002,2, (2t + 1 + cp) A K p ? l + 1 + °*>, 

Ds = a 2,2 p ( 2» + 1 + cp) A_ ,_ h _ p £" + x + <*, 

us = a 2,2,2,,2g^, q A_ i+j ^ _ p + <2 ^ 2i + CJ> , 

etc., 

where i,j, p, q take all integral values from + oo to — oo. Equate to zero the 
coefficients of the various powers of £ for all values of i, p, in order to obtain the 
requisite equations of condition between the coefficients. It is necessary that 
cp should not be an integer for any value of p ; otherwise a set of equations of 
condition would be lost and an indeterminateness would result. As far as can 
be seen, cp is not an integer for any finite values of p. With this proviso there 
results 

Z£ g [(2i+c P y-{2j+l + cq)(2i+cp—2j-cq-l)-2m(2i+cp-4j—2cq-2) ^ 

+1 wi 3 ] A 9 A )-i. q-p + 1 ™%2 q (A tAi-i-i. p-q + A ). ,^-i-i-i. _,_,)= 0, 



{2i+cp)Z£ q [2i+cp— 4j— 2cq— 2— 2m] A u q A,_ t , „_ p 

~T? m Z,jZ,q\A), q-^-i—) — \,p — q Aj, qA_ i _ i _ x< —p — q) =0, 



H3) 



where in these equations /and q take all values from + oo to — oo, and the 
equations are true for all values of i and p within the same range with a single 
exception — when i = p = , the right-hand side of the first equation is G, the 
constant of the Jacobian integral. 

It will be noticed, either from the way in which these equations are obtained 
or immediately from their form, that the coefficients in square brackets are 
directly deducible from those of Hill's equations (Vol. I, p. 133) by putting in 
the latter i + cp for i and j + cq for j. The process he gives for isolating the 
terms of principal importance in corresponding pairs of equations can therefore 
be employed here. 

Multiply the first equation by 2 and subtract the second multiplied by 3 ; 
we obtain 

2 j2« [fU> »t P, q) A j-i, i-pA i + 9«^-i-.,-i. -p-gA J = o , 
where 

f{J, i,P,q)= — (2f + cpf + 4 (2i + cp) + 2 + 4 (2t + cp + l)(2j + cq) 

+ 2 (2/ + cqf + 2m { 2* + cp + 4 (2/+ eg) + 4} + 9m 2 . 



Brown : The Elliptic Inequalities in the Lunar Theory. 251 

This single equation includes both equations by giving positive and negative 
values to i and p. The previous equations were sufficient if only positive or 
negative values were taken. A corresponding pair of equations is defined by 
the values i and p, — i and — p. The terms of principal importance for obtain- 
ing A iv and A_ it _ p are got by putting in the first of this pair of equations 
q=p,jz=i and q = ,j = 0, and in the second q=0, /=0 and q=: — p, 
j = — i. These terms are respectively then 

/(•' » *i P>P) A i,pA,o+f( 0, i, p, 0)i._(,_ p 4 , o. 
/(0, —i, —p, 0)A i>p A , +f(—i, —i, —p, —p)A_ u _ p A^ . 

Regard the equations as simultaneous ones to determine A ip and A_ i _ 1 „ 
solve so as to obtain A itP and A_ ( _ _ p separately. The common divisor is 

12 (2* + cpf [2 (2t + cp 4- 1)(2» + cp — 1) — 4m + m 2 ] , 
and either equation is included in 

2£, {(/, », jp, ?H\ «4*-«. »-? + (». i>) 4/. ?4i-i-i. P -« 

+ 

where 

(2y + cj) { ( 2t + cp) 2 — (2 4- 4m — m 2 ) + ( 2i + cp)(2 4- 2m) \ 

, • • x + (2/-+ cqf(2i + cp — 2 — 2m) 

M'V-P'^ (2t + cp)j 2 (2* + cp)» — 2 — 4m + m»"} 

r ^— 3m 2 (2i-f cp — 2)(2i + cp — 2 — 2m) — 6 —1 2m — 9m 2 

\ % ' p '~ 4(2t + cp) 2 2 (2» + cp) 3 — 2 — 4m + m 2 ' 

K 1 = 3m 2 (2 +10t + 5cp)(2t + cp— 2— 2m) + 6 4 - 12m -f- 9m 2 

L*'-P-l 4(2» + cp)* 2(2i + cp) 2 — 2 — 4/n + m 2 

In this equation the multiplier of J^^A), o or (•»> hP,p) 1S — 1 and that of 
A_ it _ p J 0i or (0, i, p, 0) is 0. By making the substitutions mentioned above, 
this equation might also have been deduced directly (Vol. I, p. 135). It is evi- 
dent that all through the transformations cp is associated with 2i and cq with 2/. 
The single equation (4) includes all the equations when i and p receive positive 
and negative values. When i = p = , the right-hand side is no longer zero, 
but depends on the constant of Jacobi's integral and connects the latter with the 
other four constants introduced into the assumed solution. 

Let us look for a moment at the way in which the primary forms are 
obtained from this general form. For the first form take the set of equations 



-j-l, P-Q \<A\ 



00 
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given by i = and p all integral values from +00 to — », and put in them 
w = 0; since they reduce to the given forms for the first solution, this entails 
j = only. In addition as n' = and g = en, n becomes equal to g and c = 1 . 
The coefficients (i, p) and [i, p] disappear, and we have for the determination 
of the Bessel's functions 

or % [> 2 + ? (p - 2) + 2p — 2] A 0t , ^ «-p = , 

(except when £> = 0, when the right-hand side is C), and the assumed solution 

becomes, when the axis of the ellipse is made to coincide with the axis of x; that 

is, when £ x = t , 

x = aoS p A 0<p cos(>+ l)g(t — 1 ), 

y = «o2,A P sin (p + l) g (t — * ) . 

Comparing, we obtain 

a o-™o, t ^o-^-o, — 2 = a i 











a oA^ 


— a A 0i 


-2 — 


:&, 


and 


as 


a first 


approximation 
















A, = 


1, A 


—i—~ 


-fe, 


A, 1 



and A 0t p is of the order e v at least. 

The particular solution is deduced by making p — and q = or e = , 
when it reduces to that given by Dr. Hill. Hence a 9 .A i>() = a i> i n hi s nota- 
tion. As a, is of the order to 2 * at least, it is evident that unless the denomina- 
tors contain powers of to as factors, A tt p is of the order m ii e p at least. 

The equations can then, provided m and e be sufficiently small, be solved by 
successive approximation, and this method is used below. In certain cases a 
power of to will appear in the denominator, but this does not affect the principle, 
as even then each approximation carries us at least two orders higher in to; 
in the solution given above it is impossible for more than the first power of to 
to be contained as a factor in the denominator. 

II. 

The system of denominators 

(2i + cjp)[2 (2t -f cpf— 2— 4m + to 2 ] 
deserves special examination. The value of c is discussed below. For the pur- 
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pose in view here, I assume that it is known and that it has the expression in 
series given by Delaunay* or the numerical value given by Dr. Hill for its prin- 
cipal part.f These are 

c= 1.071583 

c = 1 + m — 1 m* — -W- m 3 . . . . , where m = n'/(n — n') . 

Consider the algebraical expression of the two factors of this denominator- 
system when the second value of c has been substituted. The terms independent 
of m can only vanish in the first factor when p = — 2i, and in the second factor 
when p = — 2i± 1 . Hence the coefficients A ( _ 2U A ( _ 2i±1 will have their orders 
lowered by at least one power of m due to this cause. Further, they cannot have 
their orders lowered by more than one power of m from these denominators, since 
the coefficient of m cannot vanish while that independent of m vanishes for any 
values except when i=.Q, p= ± 1 , a case which does not come in here. Owing 
to the fact, however, that any coefficient A it p will depend on those of lower orders, 
this rule does not give the lowest order. It can be easily seen from the equa- 
tions (4) what the order of any coefficient is. In particular, we have 

~ . ~ . tM*. -i> -^-<.i °f order em* -1 , (i= 1, 2 . . . . <») 

— — 1 1 — "T" 1 ) 

p= + 2 i = — l|^ i( _ 2 , A_ it 2 of order eW'- 2 , (i=2, 3 ») 

p=— 2 i= -f lj^i,_ 2 , A_ hi of order e?m, 

and so on. These are well known, but the simplicity with which they are 
obtained seems worth mentioning. 

Eeturning to numerical values, the equations 

2i+ p X 1.071583 = 0, 



p = 
p 



2i+p X 1.071583 . .. . = ±\/l — 2m + m 2 /2. 

= ± 1.076303 

solved approximately in integers will give those coefficients whose expressions 
have small denominators. The set of values 

p = — 2i or — 2i ± 1 , (i = + o> . . . . — oo ) 

•Comptes Rendus, t. LXXIV, p. 19. t Acta Math., Vol. VIII. , p. 36. 
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will as before give a series of such coefficients, but the denominators increase 

with the value of i. Again 

p = ± 15, i= =F 8 

will be the smallest denominator of another series from the first equation and 

^=±14, i==F8 

from the second equation. In this way it is possible to find those terms with 
comparatively large coefficients. 

If we had included the terms due to parallax, solar eccentricity and inclina- 
tion, the denominators would have been of the form 

(i + cp -f gq + rm)\2 (i -f cp + gq -\- rmf — 2 — 4m + m i ~\ , 

where i, p, q, r receive all positive and negative integral values, and where g is 
the relation of the motion of the node to the difference of the mean motions. 
The values for which the first factor becomes small are well known, since it 
occurs directly in integrating the equations used in the older methods. The 
equation. 

i -f- cp + gq + rm = ± V 1 — 2m -f- m 2 / 2 

will have approximate solutions in integers, and those solutions in which 
i, p, q, r are not very large will determine coefficients great in comparison with 
their apparent order. These trial solutions would be obtained by inserting the 
numerical values of e, g, m. 

The terms obtained by making the first factor small are the well-known 
long-period terms ; those obtained from the second factor have periods of nearly 
a month. This closer approximation than the value unity, usually taken instead 
of the expression V 1 — 2m + »t*/2, might reveal new terms hitherto neglected. 
The monthly differences between observation and the theoretical tables are still 
sufficiently large to justify an examination of these short-period inequalities. 

III. 

Determination of the Coefficients dependent on the First Power of the Eccentricity. 

It has been shown that A Kp is at least of the order e v in e. In order 
therefore to obtain the values of the coefficients depending on the first power of 
e only, p will have the values -J- 1 and — 1 , while the corresponding values 
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which q can take are +1,0 and 0, — 1 . Put for convenience 

the equations (4) giving e ( and e'_i will become 

%j\ (J> »"i 1 » 1) «W/-i + (/ + «'f »i 1.0) efa+i ] 

+ 2(», l)^-,-! + 2 [t, 1] e;a_ i _ i _ 1 }= 0, I 

2j{ (y — », », — 1 . 0) eflj-i -f (y, — », — 1 , — 1) s-a j+i I 

+ 2 [— i, - 1] e^_ i _ 1 + 2 (— », — 1) e/a-i-,-! J = . J 

Since the quantities a ( are known, these equations will suffice to determine the 
relations of e s and sL t to e or e' . There are in all as many equations as there 
are unknowns s it sL { - In order therefore that these equations may be consistent, 
some relation must hold between the quantities a { and the coefficients ; it is 
from this relation that the value of c is to be found. For the sake of the argu- 
ment which follows, the equations will be given another form. 

In (3) put p = 1 , q = , 1 , add and subtract the results. We obtain 

2,{e*o,_ i + 0J, ««,_«+ Gj+i %i a J+i -t-iiri s (a i _j_ 1 + 5a_ i _ j _ 1 )\Xj 

+ cXj (Hj, t a»_, - HJa j+t ) Y s = 0, 

^■{c^-i + Gj^aj-i— G j+i<i a j+i -f f ^(d^.i — so..,-,.!)} r; 

+ cS, (^, i %_< + 5,'a,+ ,) -X) = , j 
where 

A , = £(2/ + l) 2 + (2/ + l)(i + 2m) + fm 2 , 

^ ii = f(2y+l) + .-+2m, fi/ = i(2y+i), 

Xj = e j + 6j, Yj — Sj Sj. 

Let oo denote the number of integers in the sequence 1, 2, 3 . . . . a>, there 
are then 4X» + 2 of these equations and the same number of unknowns 
X j and Yj. As these equations are linear, suppose the unknowns eliminated by 
a determinant, which may be looked upon as in equation in c having 4 X » -+■ 2 
roots. It is immediately evident from its form that it involves only even powers 
of c, since the coefficients G, H, H' do not contain c. 

The above equations are convergent series when |c?I^|, |clX,| are not 
infinite for any value of/, and the series 2|oi| is absolutely convergent. 

*The coefficient here called a { is therefore the same as Dr. Hill's as/ao. The change is convenient, 
since we are principally concerned with the ratios of the coefficients to a» ; <u will be a numerical quan- 
tity ultimately. 
34 
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In order to find what values of c will make these equations consistent, it is 
possible to proceed as follows. In either of equations (3) put c + 2s for c and 
i — s for i, where s is an integer, the equation considered will remain unaltered. 
Hence the result of the first substitution causes the equation with suffix i to 
become that with suffix i + s . Any such change therefore merely permutes the 
equations into one another. Hence if c be any value of c which renders the 
equations consistent, any one of the set of values 

±c , =b(c ±2), ±(c zh4), ±(c :fc6), 

will have the same effect. If we form a determinant by eliminating the 
Xj and Yj, and if this determinantal equation for c have one root c , it will 
have as roots the set of values given above. This determines 2 x <» + 1 values 
of <? when c is known. As there are 4 X °o + 2 values altogether, there are two 
values of c 2 which are both non-infinite and which do not differ by any even 
integer. It may be shown that there is one value of c\ which is zero. Dr. Hill,* 
in obtaining the other value of c 2 , submitted equations (1) to a variation h and 
obtained two equations linear in &u, 8s. A solution of these equations might 
have been obtained by putting 

and the corresponding value for 8s, that is the solution considered above when 
the second power of the eccentricity is neglected. He points out that 

8u — Du, 8s = Ds, 

where u, s represent the second particular solution is also a solution of the equa- 
tions. Now Du, Ds are odd power series in £, and therefore the above solution 
demands that c shall have an even integral value. We thus obtain the result 
that the roots of the determinant in & are 

c 2 , (Co ± 2) 2 , (c ± 4) 2 , . . . . (cb ± 2if W ,= j ; 2 . . . . oo 



2 , 2 2 , 2 2 , 4 2 , 4 2 (2if, (2if, 

The second series of values does not appear in Dr. Hill's determinant. He puts 

$u = v.Du, 8s = w.Ds 

*On the Part of the Motion of the Lunar Perigee, etc., Acta Math., Vol. VIII. 
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and forms the equations in v, w; these are reduced so that v and w can both be 
expressed by a single equation 

Z> 2 w = 0w, 

where © is a known function. The determinant which is ultimately formed 
then consists of an equation in c 2 which has only 2 x °o + 1 roots, comprising 
those of the first series above. 

It might seem that it was now possible to proceed to obtain the value of c 
as Dr. Hill has done in the memoir just referred to. When this is attempted, 
however, we are met at the outset by the difficulty that the determinant is not 
convergent. The coefficient of the highest power of c 2 in it is the square of 



= V, 



which does not satisfy the second condition given by M. Poincare, namely, that 
the sum of the absolute values of the non-diagonal elements should be finite.* 
Also it does not seem possible to find any set of factors which would render this 
determinant at once finite and adapted to calculation. 

Let A (c) denote the determinant form obtained from equations (6). We have 

(cos 7tC — l)(cos 7tC — COS 7tC ) = A'(c) -J- v 2 , 

where A' (c) has been deduced from A (c) by suitably multiplying the rows of this 
latter so as to make the coefficients of the highest power of c 2 on each side of the 
identity the same. The determinants A' (c) and v are n °t finite, and it does not 
seem possible to find any set of factors to make them both so. All that can be 
said at present is this : as c is known to have a determinate and finite value (in 
the case of our moon), the ratio A' (c)/V 2 i s no * infinite for finite values of c and 
will not be zero unless the value given to c be one of the roots of the equation 

(COS 7tC — l)(COS TtC — COS 7tC ) = . 

For the calculation of c this method appears to be of no use, but the algebraical 
forms which it brings forward are not devoid of interest. 



* Bulletin de la Societe mathematique de France, t. XIV, pp. 77-90. 
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The value of c can be calculated by continued approximation in connection 
with the coefficients e,, si when we know a first approximation to it. It is easily 
seen that the value we want is nearly unity ; from equations (5) this would give 
as a real first approximation 

c= + ^/(l-2m+^). 

As a calculation by this method would involve recomputing the coefficients in (5) at 
each step, the method would be exceedingly laborious. I have therefore thought 
it best, instead of attempting to find c, to assume Dr. Hill's value and calculate 
the corresponding values of e t and e( in terms of e and e' . There are then unused 
the two equations for e and s^, which by substitution of the values of s { and e' ( 
should give the same ratio e /e£ if the assumed value of c is correct. Since the 
equations are all linear, this is not difficult, and such a proceeding will give a 
real verification of the value of c as far as the approximations are taken. 

IV. 

Calculation of the Coefficients e 4 and e 4 '. 

In order to calculate the coefficients in the equations (5), it is convenient 
to put them into that form which will render the labor as small as possible. Let 

Li ) _ (2i±cf — (2+ 4m — m % ) + (2i±c)(2 + 2m) 
L[\~ ~~ (2i±c)\2(2i±cf — 2- 4*n + m»f ' 

M t )_ 2i±c — 2 — 2m 

M{\~ ~~ (2i±c)\2(2i±c) z — 2 — lm + m % \ ' 

where L it M ( refer to the upper signs and Li, Ml to the lower. We then have 

/(i, /, l,q)=(2j + qc) A + (2/ + qcfM it 
f(i,j, -1, q) = (2j + qc)Ll+ (2j+ qcfM{. 

When we neglect all terms except those depending on the first power of the 
eccentricity, q takes the values 1 , in the first of these expressions and the values 
0, — 1 in the second ; and the first terms in equations (5) take the forms 



where 
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Pi, i = (2/ + o) a + (ay + c) 2 if,, g,, 4 = 2jL t + 4y 2 i/;, 
p;. , = (2y- C ) a + (2y - c) 2 ^, g; , = 2/24 + 4y»jf«. 



These forms are given in extenso because the same transformations have been 
used to obtain the whole set of coefficients A p ^ { calculated. The quantities 
Li, M it L'i, Ml are calculated once for all for each equation; it is then easy to 
obtain each coefficient P, u etc. Their values are verified by the known results 
P (> * = Pii= — 1 . The coefficients once calculated, the work proceeds in the 
usual way by successive approximation, s t and si being found as linear functions 
of the unknowns s and e' . 

The coefficients a t are, as has been mentioned, the same as Dr. Hill's' a,/a , 
and the numerical values of them are taken from his paper (Vol. I, p. 247). 

It is useful to have an equation of verification at each step. Either of equa- 
tions (6) or the second of equations (2) answer this purpose. One of equations 
(6) with suffix i enables us to verify s t , sLi, e_, and s' t at once; and since the 
result of substituting the values obtained ought to render it an identity, s can 
be put equal to e« for this purpose. In this way the following table of results 
has been verified : 

1 st approx. +.01998 37455 s +.2Q559 07151 ej —.01054 25540 s —.07777 75955 e' 

2 d " +.00001 51270 s +.00008 82876 e[> —.00000 42515 e —.00001 79511 ^ 

3 d " +.00000 00038 s +. 00000 00085 4 —.00000 00003 e +.00000 00036 e^ 

£_i = +.01999 88763 e +.20567 90112 ^ s{ = —.01054 68058 e —.07779 55430^ 

<?i e-i 

1 st approx. +.00308 02608 e — .00092 87880 e^ —.00108 65604 s —.00019 44186 e^ 

2 d " +.00000 00319 e +.00000 07811 *„' —.00000 00355 e„ —.00000 15807 e^ 

3 d " +.00000 00000 s +.00000 00002 e£ —.00000 00001 6 —.00000 00006 ei 

e 1= +.00308 02927 e — .00092 80067 si ^=—.00108 65960 e — . 0001 9 59999 s' 

1 st approx. +.00001 15195 e +.00007 34781 ei —.00005 93615 s — .00043 19502 e^ 

2 d " +.00000 00010 e — .00000 00090^ —.00000 00261 s —.00000 01280 e^ 

e_ 2 = +.00001 15205 e +.00007 34691 s' s' 2 = —.00005 93876 s — .00043 20782 e^ 
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1 st approx. +.00001 47373 e — .00000 85454 s' +.00000 01043 e —.00000 08608 s' 

2 a " +.00000 00003 e +.00000 00076 e' +.00000 00000 e —.00000 00010 ^ 

e 2 = +.00001 47376 e — .00000 85378 e' sL 2 = +.00000 01043 s —.00000 08618 e^ 

e_3 = — -00000 00193 <? — .00000 01734^ £3 = — .00000 04039 e —.00000 29218 e 

e 3 = +.00000 00843 e — .00000 00708^ ei 3 = +.00000 00024 e —.00000 00055 s' 

e-i=— .00000 00001 6 — .00000 00012^ e{ = — .00000 00029 e — . 00000 00212^ 

A second approximation produces no effect in the last three coefficients. 
The advantage in point of accuracy is easily seen by the speed with which the 
coefficients approximate. 

In equations (6) put i = 0; they become 

2,- U*+ (2/ + 1 + 2m) 2 + \ m*\ a, + 1 m'a-j- J X t + oXj (2/ + 1 + 2m) a, Fj = , 
2, [c 2 a,— SmV.,-. J I) + cXj ( 2/ + 2 + 2m) a 5 X 5 - . 

Remembering that 

Xj = % + fy , Fj- = Sj Sj , 

these give on substitution of the values just obtained for f-, ej, 



whence 



2.01291 56632 7 X + 1.24333 41788 2 F = 0, 
2.31465 02085 6J. + 1.14989 92484 3 F = 0, 



^o_ 



p^ = — 2.01291 56632 7, 

y 

^ = — 2.01291 56634 5. 

The difference is probably due to accumulated errors. Taking the mean of these 

two values, 

Y 

V- = — 2.01291 56633 6. 
^0 

Hence with an error not exceeding one unit in the last place given, 

■/= —.49679 18022. 
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We may therefore say, from the way that c is involved in the equations given 
above, that the assumed value is certainly correct to the tenth place of decimals. 

As only one of the symbols e , el, X 9 , Y can be arbitrary, we must decide 
which is the most convenient. Below it is shown that if e be the value of the 
eccentricity which Delaunay has used, Y differs very little indeed from 2e. In 
the case of our moon, F will therefore be just over .1, a convenient number to 
deal with. Hence Y seems the best suited to be the arbitrary constant, and all 
the coefficients, will be expressed in terms of it. 

From the ratio X / Y we find the ratios e f Y , 4/Y , and substituting in 
the table given above, we obtain as the final values of the coefficients the 
following : 



fo 



= + .25160 40989-, 



^ = —.14889 75297, 
■«o 



Y 
si 



= _.74839 59011 



^= + .05556 82459, 

•* n 



*+l_ 



+ .00146 95307, 



£1 



-^±= — .00012 67065, 



F 

f-3 



;- = — .00005 20854, 



= +.00001 00977, 



Fo ~ 



+ .00030 84234, 



= + .00000 06713, 



,=2 = +.00000 01250, -^ = + .00000 20851 

-*0 -*0 



To 

Y 



= + .00000 00742, 
= +.00000 00009, 



-j 



= +.00000 00048, 



■1 



Let 



Y 
1=. c(n — n')(t 



+ .00000 00243. 
t ), 



D=(n — n')(t - tj , 

and let v denote the excess of the true over the mean longitude. 

The values just obtained give for the portion depending only on Y , 

r cos v = a Y [— .49679 18022 cos? 

— .09332 92838 cos ( 2D — I) + .00134 28242 cos (2D + I) 

+ .00025 63380 cos ( 4D — l)+. 00001 07690 cos (4Z> + I) 

+ .00000 22101 cos ( 6D — l) + . 00000 00790 cos (6Z) + I) 

+ .00000 00252 cos ( SD — I) + .00000 00005 cos (8Z> + I) 

+ .00000 00003 cos (lOD — l)], 
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r cos D = a F [+ 1.00000 00000 sin I 

+ .20446 57756 sin ( 2D — 7) + .00159 62372 sin (2D + I) 
+ .00036 05088 sin ( 4D — I) + .00000 94264 sin (4tD + I) 
+ .00000 19601 sin ( 6Z> — Z) + .00000 00694 sin (6D + 
+ .00000 00234 sin ( %D — I) + .00000 00004 sin (SD + I) 
+ .00000 00003 sin (10Z) — Z)]. 

The coefficients of the terms with arguments &D + I and 10D — I have been 
added by induction. 

Dr. Hill suggested (Vol. I, p. 145) that a should play the part usually 
assigned to a in the Lunar Theory. In the same way it seems best that F 
should be regarded as the arbitrary constant of this solution. The above then is 
complete so far as the first power of F is concerned. 

It will be useful to compare these expressions with the corresponding ones 
obtained by Delaunay.* For this purpose the values found above must be 
transferred to polar coordinates. We have 

tan v = y' jx'. 

In y' and x' are included the terms depending on tn only. Denote by the opera- 
tor 8 the aggregate of the new terms added, depending on the first power of F , 
then _ x'hy' — y'h x< 

x -f-y 

In this formula a/, y' denote the values of rcosu, rsinu when Y is neglected, 
and &*/, hy' the values when the terms depending only on F are taken into 
account. The value of hv is most easily found by giving to 2D the values 
0, 30, 90, 150, 180 successively and calculating the coefficients of sin I and cos I 
in the several cases. In this way are obtained the series of coefficients in longi- 
tude. That of sin I comes out to be 

.99972 87063 F . 

Delaunayf fixes the coefficient of sin I as that which would be obtained in a 
purely elliptic motion. To the first power of the eccentricity this is 2e. Hence 

2e= .99972 87063 Y , 
and therefore F = 2.00054 2735 e. 

* Memoires de l'Academie des Sciences, t. XXIX, pp. 818, 821, etc. 
t Delaunay, Mem. cit, p. 798. 
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Using Delaunay's value of e, namely, 

e—. 05489 930, 
we obtain the following series of coefficients in longitude : 

+ 4607".984 sin (2D — l) + 35".2200 sin (AD — I) 

+ 0".2906 sin (6Z> — I) + 0".0027 sin (SD— I) 
+ 174".8610 sin (2D + I) + 1".4460 sin (AD + I) 

+ 0".0121 sin (6D + I) + 0".0001 sin (8D + I). 

With the assumed value of e, these are correct to the last place of decimals given. 
The coefficient of sin (8Z> -f l) is added by induction. Delaunay's values for the 
corresponding parts are 

+ 4607.771 sin (2D — 1) + 35".1542 sin(4Z>— I) + 0".2174 sin(6Z>— I) 
+ 174".8660 sin (2D + I) + 1".4094 sin(4l> +1) + 0".0055 sin(6D + I). 

[To be continued.] 
Havekford College, December, 1892. 
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